CO depletion and deuterium fractionation in pre-stellar cores 
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ABSTRACT 

We report the detection of D2CO in a sample of starless dense cores, in which 
we previously measured the degree of CO depletion. The deuterium fractionation 
is found extremely high, [D2CO]/[H2CO] ~ 1 — 10%, similar to that reported in 
low-mass protostars. This provides convincing evidence that D2CO is formed in 
the cold pre-stellar cores, and later desorbed when the gas warms up in protostars. 
We find that the cores with the highest CO depletions have also the largest 
[D2CO]/[H2CO] ratios, supporting the theoretical prediction that deuteration 
increases with increasing CO depletion. 

Subject headings: ISM: molecules — ISM: abundance — stars: formation 

1. Introduction 

Multiply deuterated molecules have drawn much attention in the past few years, with 
the detections of D2CO (Turner 1990), ND2H (Roueff et al. 2000) and even ND3 (Lis et 
al. 2002a; van der Tak et al. 2002). D2CO was found to be particularly abundant (with 
[D2CO]/[H2CO]~ 5 — 40%) in low-mass protostellar environments (Ceccarelli et al. 1998; 
Loinard et al. 2002). The mechanisms leading to such high degrees of deuteration remain 
however debated. According to Tielens (1983) and Charnley et al. (1997), high deuterium 
fractionation can be obtained via active grain chemistry: the deuteration of the molecules, 
stuck onto the grains, is regulated by the atomic D/H ratio of the gas accreting into the 
mantles. The other major way to forming deuterated molecules is through pure gas-phase 
reactions in presence of strong CO (and other molecules) depletion (Roberts and Millar 
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2000a, hereafter RMa). In this case, the deuterium is propagated to other molecules mostly 
by reactions with H2D"'", and the [H2D"'"]/[Hjj"] ratio, which is also proportional to the atomic 
D/H ratio, plays a crucial role in determining the molecular deuterium ratio. In both ap- 
proaches, a low temperature (T < 20 K) is necessary to obtain enhanced D/H or [H2D+]/[H^] 
ratios and therefore enhanced levels of molecular deuteration. 

The gas-grain interaction mechanism was favoured to interpret the observations towards 
1RAS16293-2422 (Ceccarelli et al. 2001, and references therein). In this scheme, D atoms 
react with the CO molecules on the grain surface during the colder, pre-coUapse stage, and 
the obtained molecules arc later released into the gas phase through heating from the newly- 
formed protostars. Pure gas phase models (Roberts and Millar 2000b, hereafter RMb) could 
however well account for the observed [ND2H]/[NH3] ratios found by Roueff et al. (2000). 
Recently, D2CO was detected in L1689N, the cloud which IRAS16293 is embedded in, with 
an abundance ratio similar to that found at the border of the envelope of the protostar 
(Ceccarelli et al. 2002). This detection posed new problems, since none of the aforementioned 
models can account for the abundance of D2CO in that cloud. In order to have a high atomic 
D/H ratio, both models require strong CO and O depletion, which is not consistent with the 
high atomic O abundance observed (Caux et al. 1999; Vastel et al. 2000; Lis et al. 2001). 
Moreover, a former colder, denser (and possibly more depleted) phase cannot be invoked for 
L1689N. 

Prc-stcllar cores, as (starless) progenitors of low-mass protostars, offer the opportunity 
to study the cold, dense evolutionary stage in which the deuteration of molecules supposedly 
occurred. In addition, these objects are known to show strong degrees of CO (among other 
molecules) depletion, which is a necessary requisite for the models: the amount of deuteration 
should be hnked to the degree of CO depletion. We recently studied the degree of CO 
depletion in a sample of pre-stellar cores, and showed that the depletion degree increases 
with density (Bacmann et al. 2002). In this Letter, we report the successful search for D2CO 
in those pre-stellar cores and test quantitatively the hypothesis that stronger CO depletions 
lead to larger [D2CO]/[H2CO] ratios. We first briefly describe our multi-transition H2CO and 
D2CO observations in Section 2 and present our results on the H2CO abundance and D2CO 
over H2CO abundance ratio in Sections. We then discuss our results and their implications 
on deuteration mechanisms and finally conclude in Section 4. 

2. Observations 

We observed in August 2001 a sample of cores for which we had already discussed the 
CO depletion (Bacmann et al. 2002) in various transitions of H2CO and D2CO (cf. JPL 
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Catalog) at the IRAM 30 m telescope on Pico Veleta in Spain. The spectra were taken 
in the position switching mode at the core central position (maximum of the continuum 
emission). Additional observations (small maps around the peak) to check the extent of 
the emission of H2CO(2i2-lii) and H2CO(3i3-2i2) were carried out in November 2001. The 
rms reached were between 20 and 80 mK on the H2CO transitions and between 15 and 
50 mK on the lower energy D2CO transitions D2CO(2i2-lii) and D2CO(3i3-2i2). In order 
to estimate the opacities of the H2CO transitions, we observed in August 2002 the (2i2- 
lii) transition of H2^C0 in all cores but L1544. In all these observing runs, we used as 
backend an autocorrelator with spectral resolution of 40 kHz at 100 — 150 GHz, and 80 kHz 
at 200 — 245 GHz, yielding resolutions of ~ 0.08 — 0.11 kms~^ depending on line frequency. 
The beam and forward efficiencies at the different line frequencies can be found on the IRAM 
Web page at http://www.iram.es. Pointing was checked every 1.5 hours and found to be 
better than ^3". In the following, line fluxes are given in units of main beam brightness 
temperature. In addition to the rms noise, a calibration uncertainty of 15% at 150 GHz and 
of 20% at 230 GHz was taken into account. 



3. Results 

In 5 out of the 7 cores in our sample, D2CO was clearly detected in at least one transition 
(Table 1) ^ . The transition D2CO(2i2 — In) is a 5a — 12a detection in those 5 cores (Fig 1). 
The line fluxes given in this Table were determined by integrating the line signal between 
two given velocities (typically, where the signal is larger than the noise). For non-detections, 
the upper limit was taken as 3 times the rms on the velocity-integrated intensity, given a 
linewidth close to that of detected hues of the same molecule. 

To derive the D2CO and H2CO column densities and abundances, we used rotational 
diagrams (see Goldsmith and Langer 1999 for a full description of the method). Rotational 
diagrams give a correct estimate of the column density only when the lines are thermally 
populated and optically thin. If, however, the lines are sub-thermally excited, the rotational 
column density can be substantially different from the real one. For this reason, in this 
article wc focus only on the relative [D2CO]/[H2CO] ratio, as it is not affected by such 
uncertainty thanks to the similarity of the two molecules. Line opacity can also affect the 
rotation temperature, but has little effect on the column density ratio. This method implies 



^The non-detection of D2CO(22i — In) in L328 is questionable, since a 3a line was seen at the right 
velocity after 20 minutes of integration on both 3 mm receivers. Further integration did not confirm this 
detection. 
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identical excitation conditions in both lines. We estimated separately the ortho-H2CO and 
para-H2C0 column densities, and assumed an ortho-to-para ratio of 2 for D2CO. A single 
rotational temperature was derived simultaneously for both the ortho- and the para- H2CO 
and then used to determine the D2CO column density. In the cases where more than one 
D2CO transition was detected, the rotational temperature was derived simultaneously for 
ortho-H2CO, para-H2C0 and D2CO as well. For L1544, the rotational temperature was 
assumed to be 7.25±2.5K, which is an average of the values found for the other cores, and 
consistent with the gas temperature determined by Bacmann et al. (2002). Prom the small 
maps we made around the central position, the H2CO(3i3-2i2) emission was found to be 
extended in the cores. To correct for the difference in beam filling factors at 145 and 220 
GHz, the H2CO(3i3-2i2) was degraded to the spatial resolution of the H2CO(2i2-lii) line. 
The intensity ratio between the original and the degraded spectra (close to one) was applied 
to estimate the degraded H2CO(3o3-2o2) spectra. Note that both transitions are very close in 
frequencies and upper level energy. Finally, we corrected for the line opacity effects estimated 
by our observations of the H2'^C0 (2i2-lii) transition, which suggest opacities of 2-12 (most 
of them around 3-5). The final H2CO column densities of Table 2 have been obtained by 
taking a ^^C/^^C isotopomer abundance ratio of 60 (e.g. Boogert et al. 2000). 

Incidently, the observed H2CO line profiles are consistent (except for Oph D) with 
infalling gas, the lines being self-absorbed and the blue-shifted peak brighter than the red- 
shifted one (Leung and Brown 1977). A more detailed analysis of the observed kinematics as 
well as the absolute abundances and ortho-to-para ratios will be presented in a forthcoming 
paper. The minimum and maximum values of the column densities and the [D2CO]/[H2CO] 
abundance ratios were determined by using a Monte-Carlo method to explore the whole 
range of possible values for the rotational temperature and integrated line flux within the 
error bars. The column densities in H2CO, D2CO, H2, the D2CO abundances and the 
[D2CO]/[H2CO] abundance ratios are given in Table 2. 

D2CO is found extremely abundant in all cores where it was detected. The [D2CO]/[H2CO] 
abundance ratios given in Table 2 (between ~1% and 10%) are comparable to the values 
quoted by Loinard et al. (2002) for low-luminosity protostars. The H2CO ortho/para ratio 
was found around 2.5-3 (except for L310 where it was ~ 1.8), which is consistent with a 
gas-phase formation of H2CO (e.g. Dickens et al. 1999). 

4. Discussion &; Conclusions: Relation with CO depletion 

The high abundance of D2CO in pre-stellar cores provides convincing evidence that 
the D2CO seen in low-luminosity protostars originates in the cold, pre-stellar stage, as was 
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suggested by Ceccarelli et al. (1998) and their subsequent works (see introduction). In the 
gas-phase schemes, high deuterium fractionations of a few percent (up to ~ 40%) hke those 
measured around young low-luminosity protostars (Loinard et al. 2002) can only be achieved 
when the gas phase is highly depleted of molecules and atoms reacting with Hjj" (RMa, RMb). 
It is also the case in gas-grain chemistry models, for which high deuterium fractionations 
are obtained when the atomic D/H is high, which is in turn favoured when the gas phase 
is depleted of CO. In order to test this hypothesis, we plotted in Fig. 2 the [D2CO]/[H2CO] 
ratio against the CO depletion degree estimated in Bacmann et al. (2002) - and their Table 
2. As can be seen on the plot, the scatter and error bars are large (especially for L1544 
for which the data are not as complete), although cores with high depletion factors tend 
to have high deuteration fractions in formaldehyde. In this figure, the error bars on the 
CO depletion can be mostly attributed to the uncertainty on the dust extinction coefficient 
'^i.3mm at 1.3 mm (a factor of ~ 2). Despite the large scatter, the trend that we observe is 
consistent with an increase of deuteration with the CO depletion. To give a quantitative 
estimate, the Pearson correlation coefficient is 0.7, indicating a large probability that the 
two quantities, [D2CO]/[H2CO] and CO depletion, are indeed correlated. The models also 
state that the gas phase should be depleted of O, which is an important destroyer 
as well, only a small factor 2.5 less efficient than CO, according to the UMIST database 
(http:www.rate99.co.uk/). We have for the time being no direct information on the amount 
of atomic oxygen still in the gas phase in these cores, but some indirect indications that O is 
not much depleted indeed. On the one hand, large quantities of gaseous atomic oxygen have 
been observed in several cold, CO depleted, molecular clouds (Caux et al. 1999; Vastel et al. 
2000; Lis et al. 2002b). On the other hand, a detailed modeling of the L1544 core (also in 
our sample) molecular composition also suggests a large atomic oxygen abundance in the gas 
phase (CaseUi 2002). The apparent correlation between [D2CO]/[H2CO] and CO depletion 
and the very high fractionation observed may suggest either a similar degree of oxygen 
depletion with respect to CO (which would disagree with the aforementioned arguments) or 
a lower role of O in the fractionation (yet not supported by theoretical considerations). In 
the last case, the relatively large scatter observed in Fig. 2 may be due to the presence of O 
in the gas phase. 

The question remains whether D2CO is formed in the gas phase (Roberts and Millar 
2000b) and frozen out onto the grains or directly formed on the grains as suggested by Tielens 
(1983). If D2CO is formed on the grains, there must exist an efficient desorption mechanism 
to release part of the D2CO into the gas phase, even in such cold (T ~ 10 K) environments 
as pre-stellar cores. The depletion timescale (~ 10000 yr - to be compared with estimated 
lifetimes of 10^ yr, see Beichman et al. 1986) derived in Bacmann et al. (2002) suggests 
that desorption mechanisms indeed exist, since it is highly unlikely that all cores from the 
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sample have only reached 1% of their hfetime. The gas-grain model is thus not ruled out by 
observations of D2CO in the gas phase, even in objects with T ~ 10 K. For the rest, the same 
remarks as on the cloud L1689N (Ceccarelli et al. 2002) apply here. In the absence of O 
depletion, the atomic D/H ratio is predicted to be around 0.01 (RMa), which is not sufficient 
to account for the high levels of deuteration observed here (to account for a [D2CO]/[H2CO] 
ratio of 0.05, the atomic D/H ratio would have to be as high as 0.3 - Caselli et al. 2002b). The 
dashed curve in Fig. 2 represents the predictions of the gas-phase RMb model, which assumes 
relatively large O depletion as well. Our data compare rather well with those predictions, 
with the exception of Oph D, which shows a larger [D2CO]/[H2CO] (~ 10%) than predicted. 
However, we point out that another key parameter in the Roberts & Millar model is the 
electron fractional abundance. The predictions reported in Fig. 2 have been taken from the 
Fig. 3 of RMb, who assumed an electron fractional abundance around 10~^. This is much 
larger than the value of X{e) ~ 10"^ actually found by Caselli et al. (2002a) in L1544. Fig. 5 
of RMa suggests that larger [D2CO]/[H2CO] ratios could be reached if the ionisation fraction 
X{e) is lowered. It is also worth mentioning that the absolute D2CO abundances we measured 
for the degrees of CO depletion found in the cores (Table 2) are generally smaller than those 
calculated by RMb (between 10"^'^ for low depletion factors and 1.6x10"^^ for depletion 
factors of 16), though in the end, we find deuterium fractionations similar or higher than 
those predicted by the model. This discrepancy is strongly marked (one order of magnitude) 
for the cores with little depletion (L1689B, L1709A and L328). To summarize, gas-phase 
models, though they can account for the relative deuterium fractionation in most of our 
objects, cannot explain all of our data. Key points that still need to be fully understood are 
the depletion of atomic oxygen and its role, as well as the ionization fraction in pre-stellar 
cores. 

A.B. acknowledges support by the Cerman Deutsche Forschungsgemeinschaft, DFG in 
the course of this work. We thank the IRAM staff for their help during the observations and 
the Time Allocation Committee for awarding of the telescope time. 
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Fig. 1.— H2CO(2i2 - 111), Hi3CO(2i2 - In) and D2CO(2i2 - In) lines in the pre-stellar 
cores of the sample. The D2CO line is detected with a good signal-to-noise ratio in L1544, 
L1689B, L1709A, L429 and OphD. H2CO is double-peaked in 6 of the cores. The /sr velocity 
is shown as dashed lines. 
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Fig. 2. — Comparison of [D2CO]/[H2CO] abundance ratio and CO depletion factor (from 
Bacmann et al. (2002) for the dense cores in the sample. The [D2CO]/[H2CO] vs. CO 
depletion relation derived from the model of RMb is represented by the dashed curve. 
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Table 1. Integrated areas of lines observed (K km s ^) for the 7 cores of the sample. 
Upper limits are given for 3cr (see text). A dash means that no observation was performed. 



Core 




L1544 


L1689B 


L1709A 


L310 


L328 


L429 


Oph D 


H2CO(2i2 - 


111) 


0.62 ±0.11 


1.63 ±0.28 


1.24 ±0.22 


0.40 ±0.08 


1.52 ±0.25 


0.93 ±0.17 


1.43 ± 0.24 


H2CO(2o2 - 


loi) 




1.69 ±0.30 


1.16 ± 0.20 


0.45 ± 0.11 


1.23 ±0.22 


0.66 ±0.13 


1.40 ±0.24 


H2CO(3o3 - 


202) 


0.26 ±0.07 


1.04 ± 0.28 


0.59 ±0.14 


0.31 ±0.12 


0.34 ±0.10 


0.31 ± 0.08 


0.58 ±0.13 


H2CO(322 - 


221 ) 


< 0.083 


< 0.18 


< 0.068 


< 0.237 


< 0.090 


< 0.054 


< 0.076 


H2CO(3i3 - 


2l2) 




1.46 ±0.37 


1.10 ±0.27 


0.30 ± 0.09 


0.86 ±0.21 


0.61 ±0.16 


0.96 ± 0.23 


D2CO(2i2 - 


111) 


0.052 ±0.012 


0.058 ±0.012 


0.11 ±0.02 


< 0.046 


< 0.020 


0.071 ±0.012 


0.12 ± 0.02 


D2CO(3i3 - 


2l2) 




< 0.14 


0.089 ± 0.032 


< 0.18 




0.029 ±0.013 


0.14 ±0.04 


D2CO(4i3 - 


3l2) 




< 0.24 




< 0.16 


< 0.25 


< 0.061 


0.069 ± 0.033 


D2CO(4o4 - 


303) 


< 0.05 


0.11 ±0.05 


< 0.094 


< 0.27 


< 0.065 


< 0.043 


< 0.22 
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0.33 ± 0.02 


0.13 ±0.02 


< 0.04 


0.10 ±0.01 


0.07 ±0.01 


0.04 ± 0.01 
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Table 2. Rotational temperatures, H2CO and D2CO column densities observed in a 17" 

beam (beamsize at the frequency of the H2CO(2i2 — In) transition), as well as 
[D2CO]/[H2CO] abundance ratios, H2 column densities (from Bacmann et al. 2002) and 
D2CO abundances Xoaco (-^^0200 — -^D2Co/-^H2)- The H2CO column density values are 

corrected for opacity. Upper limits are 3a. Figures follow the convention a{b) — a x 10^. 
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-11) 


L328 




4.3li1 (13) 
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9.5 (22) 


< 6.1( 
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1.4 (23) 


1.7(- 
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